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Abstract 

We consider the soliton solutions in 1- and (l+l)-dimensional Toda lattice models 
with a boundary. We make use of the solutions already known on a full line by means 
of the Hirota's method. We explicitly construct the solutions satisfying the boundary 
conditions. The Zoo-symmetric boundary condition can be introduced by the two-soliton 
solutions naturally. 
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1 Introduction 



Integrable models on a half line have been widely studied recently. To investigate the inte- 
grable models on a half line, two approaches are already known. One is the algebraic approach 
by the boundary Yang-Baxter formalism[|l| and the other is the field theoretical one with the 
Lagrangian accompanied by the proper boundary terms 0]. The model which is integrable 
on a full line needs the proper boundary terms in order to be integrable also on a half line. 
Therefore, the boundary terms should be determined by the requirement that the infinite 
number of integrals of motion on a full line must survive after introducing the boundary. 
Both approaches assume that the particle content and the S-matrix are the same as those 
on a full line. Though they offer quite beautiful descriptions of the integrable models on a 
half line, especially in sine-Gordon model || |[] and afline Toda field theory (ATFT) ||, the 
relations between the two approaches are to be clarified. As for the sine-Gordon model, the 
semi-classical analysis has already been applied to determine the relation of the boundary 
term and the boundary S-matrix ||. 

In the Lagrangian approach, there are two ways to demand the infinite number of integrals 
of motion to survive after introducing the boundary. In the first, the given boundary terms 
determine the boundary condition of the variables or fields through the Euler-Lagrange equa- 
tion. These boundary terms give the boundary condition of the conserved currents, which 
have the same forms as those in a full line. Therefore it is possible to define the integrals 
of the conserved currents over the half line, which do not have the asymptotic condition. If 
we demand those integrals to be time-independent, we can determine the boundary terms 
inversely. The second way is to use the transfer matrix. We make the proper image of the 
half line at the boundary by adding the proper boundary terms to the Lagrangian, which 
enable us to consider the trace over the spatial direction. Therefore we can define the transfer 
matrix on a half line properly and demand it to be time- independent. 

On the other hand, Hirota's method has been used to obtain the soliton solutions of the 
equation of motion in a variety of integrable models on a full line in 1- or (l+l)-dimensions, for 
example, the nonlinear Schrodinger model, sine-Gordon model, KdV model and Toda lattice 
model0. This method gives the solutions in quite an explicit form so it is an appropriate 
method to consider the boundary conditions. 

In this paper, we discuss the classical solutions of the equation of motion with the in- 
tegrable boundary conditions in 1-dimensional Toda lattice model and the extension to the 
(l+l)-dimensional field theory (Toda lattice field theory). As for the other integrable models 
on a half line, the classical solutions are already known by the classical inverse scattering 
method^, Therefore our analysis of Toda lattice models is considered as a first step. 

This paper is organized as follows. In section 2 the Toda lattice (TL) model is investigated. 
The TL model is considered as the multi-particle systems or the model of an electric circuit 
with solitons. The equation of motion has already been solved by means of Hirota's method 
J7J. On the other hand, the boundary terms to make the TL model integrable on a half line 
have been determined by means of the transfer matrix method0. We will consider the soliton 
solutions obtained by Hirota's method on a full line and satisfying the boundary condition 
explicitly. We will see the meaning of the parameter in the boundary term in several simple 
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examples. In section 3 the (l+l)-dimensional Toda lattice field theory (TLFT) is introduced 
to clarify the notations. The TLFT is an extension of the TL model to a field theory in (1+1)- 
dimension and has the soliton solutions unlike the 0$ ATFT with a real coupling constant. In 
section 4 the TLFT with a proper boundary term on a half line is investigated. There seem to 
be several ways to introduce the boundary in this model, for example, the half number of the 
fields should vanish. However, we consider the way to introduce the boundary in the spatial 
direction, which has the connection with the (% ATFT analyzed already, therefore it is the 
most suitable to consider. We will see that the boundary terms conserving the Zoo-symmetry 
can be induced by the two-soliton solutions. In section 5 we present a summary of this paper. 



2 Integrable Boundary Condition in Toda Lattice Model 

It is well-known that the one-dimensional TL model is integrable, which can be shown for 
example in terms of the construction by the Lax pair. If we consider the integrable models 
on a half line, integrability imposes the restricting condition on the boundary terms, which 
determines the boundary condition of the model. The integrable boundary condition (IBC), 
which is compatible with the integrability, can be determined by the requirement that the 
integral of motion on a full line must survive also on a half line. In turn, this can be realized 
for example by imposing that the transfer matrix, which is validly defined on a half line, to 
be time-independent like that without the boundary. The IBC in TL model is determined in 
Sklyanin's paper0 in this scheme. In this section, we consider the classical solutions with 
the IBC according to Hirota's method. 

To begin with, we consider the TL model on a full line. The one-dimensional TL model 
is defined with the canonical coordinates (7r n ,g n ), which are functions of the time variable t 
and on a full line, integer n G (— oo, oo). The hamiltonian and the Poisson brackets are 

oo 2 oo 

H = E Y+ E e9 ' l+1 " 9n > K»,?n} = <W (i) 

n=— oo n=—oo 

Now we write down the equation of motion in terms of new variables r n = q n +i — q n . The 
equation of motion is 

r n = e rn+1 -2e r " + e r — 1 , n = 0, ±1, ±2, ±3, . . . , (2) 

here ' represents the time derivative 4. We can obtain the soliton solutions of the equation 
of motion (^) by means of Hirota's method. To apply Hirota's method, we introduce the new 
variables V n and F n as r n = ln(l + V n ) and V n = 4$ lnF n . With these variables V n and F n , 
the equation of motion @ can be written in a simpler form 

F n F n — F% = F n+ iF n _x — F%. (3) 

The above equation of motion (||) has s-soliton solutions 

F n = E eX P E faVi + E A i0^3 > 
jit=0,l \i=l i<j J 

rji = pin - Vtit + 5i, (4) 
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where <Vs are arbitrary real parameters and //j (z = 1, 2, . . . , s) is or 1. The real-valued 
momenta and the coefficient must satisfy the condition 

Vti = 2ejsinhy, e* = ±1, 



{ ™&j{Pi-Pj) \ 2 e . e .> 
Y sinh | (pi -\- Pj ) J 1 3 " 

f cosh Kpj-p^ Y 
\ v cosh i (p i +Pj)J 



(5) 



The above soliton solutions approach zero as |n| — > oo and does not diverge anywhere provided 
Pi, Cj and 5j are real numbers [[□], [T^|. 

We are ready to consider the TL model with the IBC. The hamiltonian takes the form 

oo 2 oo f a \ 

# = X! if + £ e 9 " +1 - g " + «e 91 + §e 2?1 . (6) 



a 2 



n=l 



It can be shown that the model given by the above hamiltonian is also integrable with the 
arbitrary parameters a and (3 because it has an infinite number of the integrals of motion|0. 
The equation of motion can be written down with r n = q n +\ — q n (n = 1,2, . . .) and q\ as 

r n = e r " +1 - 2e r " + e r - 1 n = 2, 3, • • • , (7) 
fx = e r2 - 2e n + (ae qi + (3e 2qi ). (8) 

Eqs.(^) for r n (n = 2, 3, . . .) are the same as those on a full line (0) and Eq.(H) for r x includes 
the variable q±. We make a remark that if lim^oo | Y%=i r i\ < 00 ? which is true for the soliton 
solution (§,?! = (qx - &) + (ft - ?a) + • • • = - ££i 1- Comparing Eqs.(@) and (0)(|), we 
obtain the boundary condition in Eq.(H), which enables the soliton solution (Q) on a full line 
to satisfy Eqs. (|7]) and (|S|) on a half line in the form 

oo / oo \ 2 

e r ° = a J] e" r * + /3 JJ e~ n . (9) 

i=l \i=l / 

In terms of F n , this can be written down simply as 

F-! = aF + /3F V (10) 

Firstly, we will concentrate the special case with a + (3 = l,(3>0 for simplicity. In this case, 
we can find two-soliton solutions with arbitrary momentum satisfying the boundary condition 
fllCf) . We consider the two-soliton solution 

-Q.t(jpn+5i , „-pn+5 2 \ , f nnch P_\ ^ _-2nt+5i+<5 2 



F n = 1 + e ~ s "(e pn+ 1 + e~ pn+d2 ) + cosh 



= 2 sinh |, (11) 
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with an arbitrary momentum p. The boundary condition (|i~0|) gives the equation of the initial 
displacements 8\ and 8 2 

- 1±p±1 (12) 



+ 13 



which is positive-definite even if a = 1 — /3 is negative. Therefore the system is stable JTO 
This equation means the following. We assume Q > 0,p > 0. In the limit t — > — oo, there 
exist two solitons with momenta (£l,p) and (fl, —p). As the time t grows to oo, F n goes to 
1. Therefore if the initial displacements of the two solitons, 8\ and 5 2 , satisfy Eq. (|T2|) , these 
two solitons undergo a pair-annihilation due to the boundary. In particular, if a = and 
(3 = 1, Eq.(|T2D is simplified as 8\ = 5 2 . It means that in this case, the shift of the initial 
displacements due to the boundary is independent of the momentum p of the soliton. 

Before closing this section, we briefly mention the general case with arbitrary a and (3. If 
a + (3 ^ 1, it is much more difficult to find the solutions satisfying the boundary condition 
fllQf) than those with a + (3 = 1. We consider the case in which the initial momentum p takes 
a special value so that the two-soliton solutions may exist. We consider a two-soliton solution 



-2 



-2pn+8i+S2 



= 2sinh|. (13) 

If we assume the region of (3 as — e~ p < (3 < —e~~ 2p for p > or — e~ 2p < /3 < — e~ p for p < 0, 
this soliton-solution satisfies the boundary condition (|1(J) only if 

e- p -a + \(3\e p = 0, (14) 

S 1+ S 2 _ Iff I ~ e " P n 

where the value (p!5|) is made positive. We assume Q > 0,p > as before. Considering the 
limits t — > ±oo, the initial soliton with momentum (fi,p) is scattered by the boundary to the 
soliton with momentum (— 0,p) if the momentum p satisfies Eq.(|T4D and the initial phase S 2 
are determined from 8% by Eq . (|15D . 



3 Definitions and Solutions in (l + l)-dimensional Toda 
Lattice Field Theory 

In this section, we consider the extension of one-dimensional TL model to a (1+1)- dimensional 
field theory (TLFT). As well as other (l+l)-dimensional integrable classical field theories, the 
equation of motion of the TLFT can also be solved by means of Hirota's method. (1+1)- 
dimensional TLFT (a^ Afflne Toda field theory with the real coupling constant) is defined 
by the action 

/CO (-OO (-OO poo I 1 _, 1 _, °° ^\ 
da / drC = / da dr -(«9 T 0) 2 - -(«9 CT 0) 2 - ^ e 5 ^ . (16) 
-oo J —oo J —oo J —oo \ 2 Z . / 
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Here we explain the notations briefly. Let R°° be the infinite-dimensional vector space. 
0(<t, r) G R°° is an infinite-dimensional vector field and c?j = ej + i — e*i, where {c{\ is the 
standard basis of R°°, is a simple root of Lie algebra. If we put di • (f) = (f>i, the equation 
of motion can be written as 

dl<j) n -d 2 a (j) n = 2e*" -e " +1 -e*"" 1 , n = 0, ±1, ±2, • • • . (17) 

We carry out the calculation in the light-cone coordinate z = (a — r)/2 and z = (er + r)/2 
(9 = <9 2 ,<9 = (9 2 ) to keep the notation simple. If we introduce new variables F n and V n 



4> n = ln(l + V n ) and = — dd\nF n as before, the equation of motion of the model (|T6| ) can 
be written in terms of F n 's as 

dF n dF n - F n ddF n = F n+1 F n _! - F 2 . (18) 

The equation of motion (flSf) has the s-soliton solutions 

= Z exp ^ /iiPi + AijUiHj 

£i=0,l \i=l i<j 

Pi = cijZ + + cin + 8i, (19) 

here the initial phases 5i are arbitrary and the momenta cii,bi and q and the coefficients Ay- 
must satisfy the condition 

a-ibi = — 4sinh 2 — , 
2 



a j)(h — bj) + 4 sinh 



2 Ci Cj 



e 3 = cTF" ( 2 °) 

(ai + ai )(6i + 6j) + 4 sinh 2 - J 



Several remarks are in order. Firstly, from the requirement that the solutions given by ([jj 
do not diverge on a full line, i.e. V n > — 1, we will get the condition that e AiJ > and that 
Si's are real values. Secondly, the above given solutions are soliton solutions. It is well-known 
that contrary to the sine-Gordon case, there is no soliton solutions in the (real-coupling) c% 
ATFT for finite N, for example sinh-Gordon model. The reason is that in the case of a$ 
ATFT for the finite N, the cyclic condition in the direction of the affine coordinate, that is 
i + N = i, demands the momenta q to take certain complex values q = (2iri/N)Z. This fact 



causes the naive soliton solutions given in Eq.(19) to diverge to infinity in the limit a — > ±oo, 
therefore these models do not have soliton solutions. The TLFT, in spite of corresponding to 
the ATFT with the real coupling constant, has no cyclic condition in the affine direction 
so we can choose Cj to be real values, which do not cause the divergence. Conversely, we 
cannot obtain the solutions in TLFT by merely taking the limit of iV — > oo in off ATFT. 
We need to take the analytic continuation in the affine direction. 
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4 Toda Lattice Field Theory with Boundary 



In this section, we will consider the TLFT on a half line. We choose the space coordinate a 
in (—00, 0) and the time coordinate r in (—00, 00). The action of the TLFT on a half line is 
given by 

S = da drC+ dr£ B , (21) 



here the bulk Lagrangian density £ is the same as that given before on a full line (|16|) and 
the boundary term £5 is given as 

00 

C B = - E Ae^U, (22) 

i=— 00 

with the parameters A4. It is well-known that the a^-ATFT with a finite N > 2 defined by 
the above action ( |2"2] ) is integrable if |^4*| = or 2 for all i0. In these cases, it is shown that 
the integrals of motion on a full line can survive after introducing the boundary at a = 0. 
Now we consider the classical (soliton) solutions in this model on a half line. The equation 
of motion is the same as that on a full line ( |i~7|) and the boundary term £q in the action (pi ) 
induces the boundary condition of <fi at a = 



d„ 



(7=0 



2A n e K/2 - A n+1 e^ +l/2 - A-ie - l/2 . (23) 



For the later convenience, we write down the boundary condition (|23| ) in terms of the variable 
F n , which is defined the same as that on a full line, 

LK _ F'n+i _ 1 _ A ^/F n +iF n -i _ A n +i V F n+2 F n _ A n -i VF w _ 2 F n \ 1 

\ -^n ^n+i F n _i ) a ~° \ n F n 2 F n+ i 2 F n _! j e7 ~ 0, 

(24) 

here ' means the spatial derivative 8/ da. We will consider the soliton solutions (|19D obtained 
on a full line, some of which satisfy the boundary condition fl2~4]). 

Firstly, we consider the simplest example; Ai = for all i. If a solution F n is an even 



function in a, it satisfies the boundary condition (24). In particular, if we consider two-soliton 
solutions (s = 2) in (1191) with 



pi = pa — er + cn + 6, p 2 = — V a — er + era + 5, (25) 
the coefficient Ayi takes quite a simple form as 



e-^l-^smh^-JJ , (26) 

from the condition (0). Providing this solution does not diverge on a half line, the coefficient 
) must be positive, that is \p\ > 2sinh 2 |. 

Secondly, we consider the general boundary condition given by (|23|). We assume two- 




soliton solutions (s = 2) given by ([19]). Because the lefthand-side of Eq.(gJ) is a rational 
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function of e Pl and e p2 , we must impose that the square root V^+i7^iU=o takes the form 
of a certain polynomial in e pi and e P2 . This condition makes the calculation much simpler 
and gives the constraints of the momenta as 

p 2 = pa + er - cn + 8 2 , 

(27) 




where momenta p, e and c are arbitrary real values with the condition ( |2"0"| ) and the phase 
displacements 8\ and <5 2 must satisfy the equation 

l + 7 (fsmh 2 f) 

Due to the condition of the absence of the divergence, i.e. V n > —1, we must insist \p\ > 
2 sinh 2 | as before. With the above preparation, we can consider the boundary condition (p4|). 
Because we can easily obtain 



F' n + j^F n+1 F n ^ - (p - 7 )FnJ | CT=0 = 0, (29) 

we see that .4, = 27 for all z. Therefore, if we consider the two-soliton solutions, we automat- 
ically obtain the Zoo-preserving boundary condition, A4 = 2 or —2 for all i. 

We note that if only the soliton solutions are considered, it would be quite difficult to find 
the solutions satisfying the Zoo-symmetry breaking boundary condition, for example Aq = — 2 
and Ai = 2 for i 7^ 0. In particular, even if we introduce the three-soliton solutions with 

pi = pa - er + cn + 81, p 2 = pa + er - cn + 5 2 , 

p 3 = ±2a + inn + 5 3 , (30) 



which is the unique choice for the three-soliton solution to satisfy the boundary condition(p3|) 
and to have no divergence, we obtain the same result as that in the case of the two-soliton 
solutions (P9"|). It means that the third soliton given by (^) has no effect on the boundary 
condition. 



5 Summary 

In this paper, we have considered the TL model and TLFT with boundary and seen the 
meaning of the parameters in the boundary term. A few comments are in order. Firstly, 
the analysis in this paper can be applied only to the other models with the soliton solutions. 
Therefore, it is impossible to consider the ajy ATFT with finite iV and real coupling constant, 
which do not have the soliton solutions in the same way. However, it is possible to consider 
the ATFT with imaginary coupling constant. Especially, the analysis of TLFT with the 
Neumann-type boundary condition can be similarly used in that for 0$ ATFT with imaginary 
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coupling constant. Secondly, we did not analyze TL with arbitrary a and (3 nor TLFT with 
general AiS. To analyze these models, we must take multi-soliton solutions into account. 
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